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Abstract. Strong gradients in plasma flows play a major role 
in space and astrophysical plasmas. A typical situation is that 
a static plasma equilibrium is surrounded by a plasma flow, 
which can lead to strong plasma flow gradients at the separa- 
trices between field lines with different magnetic topologies, 
e.g., planetary magnetospheres, helmet streamers in the solar 
corona, or at the boundary between the heliosphere and inter- 
stellar medium. Within this work we make a first step to un- 
derstand the influence of these flows towards the occurrence 
of current sheets in a stationary state situation. We concen- 
trate here on incompressible plasma flows and 2D equilibria, 
which allow us to find analytic solutions of the stationary 
magnetohydrodynamics equations (SMHD). First we solve 
the magnetohydrostatic (MHS) equations with the help of a 
Grad-Shafranov equation and then we transform these static 
equilibria into a stationary state with plasma flow. We are in 
particular interested to study SMHD-equilibria with strong 
plasma flow gradients perpendicular to separatrices. We find 
that induced thin current sheets occur naturally in such situ- 
ations. The strength of the induced currents depend on the 
Alfven Mach number and its gradient, and on the magnetic 
field. 



1 Introduction 

Plasma flows around separatrices play an important role in 
many astrophysical and space plasmas. Significant flows oc- 
cur thereby mainly on open field lines, while the plasma on 
closed field lines is approximately at rest. One example is 
the magnetosphere surrounded by the solar wind flow, where 
both regions are separated by the magnetopause. Another ex- 
ample are coronal helmet streamers, where the closed arcade 



type magnetic structures are surrounded by open magnetic 
field lines on which the stationary solar wind is streaming. 

A third example, but with different constraints, is the re- 
gion far away from a star, which is embedded in the coun- 
terflowing interstellar medium. Due to the interaction be- 
tween the stellar wind, e.g., the solar wind, and the counter- 
flowing interstellar medium, a separatrix forms, separating 
the 'inner' stellar wind from the 'outer' interstellar medium: 
the domain inside is called an astrosphere (heliosphere for 
the sun), the corresponding separatrix is called an astropause 
(heliopause). This is a similar situation as described for the 
helmet streamers in the paragraph before, namely different 
flow regimes inside and outside of some boundary layer ('in- 
ner' and 'outer' field lines). In the case of astrospheres/of the 
heliosphere, however, the scenario implies a structure with 
almost completely open field lines. Additionally, the flow is 
non-zero also on the 'inner' field lines, but shows a strong 



gradient (see, e.g., jBaranov et aL 1970 Baranov and Kras 



|nobaev| [TWT] [Nickeler et al.||2006| for details). 

These situations where regions with and without plasma 
flow are separated by rather thin boundary layers necessarily 
lead to strong flow gradients in these layers. Within this work 
we aim to study the relation of these flow gradients to cur- 
rent sheets. Thin current sheets are important, because due 
to current driven micro-instabilities a fine resistivity occurs 
in these regions and the usual assumption of an ideal con- 
ducting space plasma breaks down. Consequently resistive 
plasma instabilities like magnetic reconnection can occur as a 
consequence of current sheets and when additional sufficient 
free energy for eruptions is available in the configuration (see 
for example Schindler and Birn 1993 1. Such processes oc- 
cur, e.g., as substorms in magnetospheres and flares in the 

for a recent study 



solar corona, (see Birn and Hesse 2009 
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of similarities and difference between substorms and flares). 
IWiegelmann and Schindler] ( [T995] l, and [Becker etaL] ( |200T] i 
studied the formation of thin current sheets as a sequence of 
quasi-static magnetotail equilibria. Quasi-static means that 



2 



Nickeler et al.: MHD with plasma flow 



dynamical effects and the influence of plasma flow can be ne- 
glected compared to magnetic forces and pressure gradients. 
This assumption is well fulfilled in the magnetosphere during 
quiet times ( jSchindler and Bim]|1982[ ) and the strongest cur- 
rents form in the center of the magnetospheric plasma sheet. 
The influence of a field-aligned parallel plasma flow in mag- 
netospheres has been studied for example in Birn ( 1991[l and 



for solar MHD-equilib ria in|Wie gelmann et al. ( 1998[), 



Petrie 



land Neukirch] ( fT999l ), [Petrie et al., (2002\ and lPetrie et al 



( 20051 1. ^ significant influence of a smooth plasma flow itself 
onto the magnetic field and plasma configuration is rather 
low for flow speed well below the Alfven speed u <C wa- Far 
less studied has been the influence of significant small scale 
gradients in the plasma flow, which is the topic of this paper. 
We are in particular interested to investigate to which extend 
thin current sheets at boundary layers can be associated with 
corresponding gradients in the plasma flow. We outline the 
paper as follows. In section |2] we present the basic equa- 
tions, outline how static equilibria can be transformed into 
stationary ones, and discuss the relation between sharp flow 
gradients and the occurrence of thin current sheets. Sections 
|3] and |4] contain applications to solar coronal and magneto- 
spheric structures, respectively. Finally we summarize our 
results in section |5] 



2 Basic assumptions and equations 

We apply the theory of ideal MHD. This is justified, be- 
cause of the high conductivity in many space and astrophys- 
ical plasmas. For simplicity and to concentrate on essen- 
tial 'flow' effects, we restrict our research to configurations 
with an incompressible plasma flow. We are interested in the 
physical effects that occur for large gradients of the Alfven 
Mach number perpendicular to the field lines. We are aware 
that compressible effects might become important for some 
space plasma applications, in particular for stratified plasmas 
like the solar corona and chromosphere, especially on large 
scales, see, e.g. [Petrie et al.| ( |2002| i or [Petrie et aL[ ( [2005[ l. 

For stationary, ideal and incompressible MHD we have 
to solve the following equations: mass continuity equation 
([T]i, the Euler or momentum equation with isotropic pressure 
P (|2]i, the stationary induction equation including the ideal 
Ohm's law ([3]), Ampere's law Q, the solenoidal condition 
for the magnetic field (jSj, and the condition for incompress- 
ibility (j6|l: 
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where p is the mass density, v the plasma velocity, B the 
magnetic field, j the current density vector, and P the ther- 
mal or plasma pressure. 

Due to the incompressibility, the mass continuity equation 
can be written as v • Vp = 0, so that the density is constant 
on streamlines. We now introduce the auxilliary flow vector 
or streaming vector w := ^pv. With the Bernoulli pressure. 



defined by 11 : 
tions as 



P 



we can rewrite the above equa- 



V-w = 0, 



1 

Mo 



{V xB)xB-{V xw)xw = Vn, 
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(9) 
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The momentum equation Eq. ([8]l is written such that the anal- 
ogy with magnetohydrostatic equilibria (MHS), given by 



Mo 



-^{V xB)xB-- 



VP, 



(11) 



is obvious. 

The assumption of a field-aligned flow enhances the prob- 
ability that the flow is stable in the frame of ideal MHD, 
(see the discussion in Hameiri' '1998). The assumption 
of V X B — leads to a vanishing electric field in ideal 
MHD. This can be seen with the help of the uncurled in- 
duction equation, Eq. (|3]l, which is basically ideal Ohm's law 
E + w X B = (E =electric field). Therefore the station- 
ary ideal Ohm's law (j9]) is fulfilled identically. Under these 
assumptions the set of equations reduces to: 



B ■ VMa ^ 

(1- 



vn = 

V-B = 0, 



= 0, 

- MVj {V X B) X B 



2^0 



V(l-Mi 



where Ma is the Alfven Mach number, defined via 

w = ±MaB/^, 



(12) 

(13) 
(14) 

(15) 



and where the ± indicates that if the pair {w,B) is a solu- 
tion of the Eqs. ^ and ( 13 i, then {—w,B) is also a solution, 
and basically also {w,—B) and {~w,—B). We will mainly 
parameterize this behaviour by Ma and by B, starting from 
the viewpoint of MHS theory, as the form ( [T3] l of the in- 
compressible SMHD equations allows to derive transforma- 
tion equations which transform MHS equilibria into station- 
ary ones (see ,Gebhardt and Kiessling 1992| Nickeler et al. 



[2006[ for the mathematical details.). 
2.1 2D stationary states 

The theory explained so far is general. In the following we 
concentrate on configurations with one axis of symmetry, e.g. 
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the z-axis in a cartesian coordinate system x,y,z. Conse- 
quently all quantities are functions of x and y only and we 
can solve the solenoidal equation ( 14i by introducing a flux 

This reduces the 



function a{x,y) of the form B — Va x e 
stationary incompressible equations in 2D to 



Ma 



vn = 



MA{a) 
(1- 



-Af|) AaVa 



\Va\_ 
2^0 



(16) 

(17) 



Now we perform a transformation by assuming that a{x,y) 
is a function of another 'flux function' A(a;,?;), i.e. a = a{A), 
such that the stationary equation, Eq. ([TTji, reduces to a 
form of the equation mathematically similar to the Grad- 
Shafranov equation 



dPMHS 

dA 



-AA. 



Mo 



(18) 



describing MHS equilibria. 

Then the equivalence between the MHS equation and the 
momentum equation of ideal, stationary but non-static MHD 
with incompressible, field-aligned flow is shown. As a is a 
function of A, Ma is also a function of the 'new' flux func- 
tion A. Derivatives with respect to A will now and in the 
following be expressed by a prime, e.g., da/dA — a' . With 
the help of the relation Ma = Ma (A) we can rewrite the Eu- 
ler equation Eq. ( [T7] i 



vn 



1 - Ml) ( a" {VAf + a'AA) a'VA 



2/10 



Mo 
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Mo 
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ZfJ.0 



(19) 



Let us remark that the Alfven Mach number Ma can be ex- 
pressed as a function of a or A, but is not restricted further. 
We have therefore the freedom to choose this function ar- 
bitrary without loss of generality. A reasonable choice to 
ehminate the term ( VA)-^ in ( 19 1 is 



(1- 



1, 



(20) 



and therefore the Euler equation ( 19 1 simplifies to a single 
partial differential equation for the new flux function A{x,y) 



vn = 



1 

Mo 



-AAVA 



dA 
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Mo 



-AA. 



(21) 



In any case, Eq. (21i is mathematically identical with 
Eq. ( frs) , but contains plasma flow. Physically this equation 
reduces to the static Grad-Shafranov equation only for the 
limit Ma 0, implying n(A) P{A). Consequently, any 



solution A{x^y) of the 'MHS' equation (21 1 (or equivalently 



the MHS equation ( 18 i) can be used to derive a solution of 
the stationary, incompressible MHD by integrating equation 
([20ll 



a = ± 



dA 



^1~Ma{AY 



(22) 



With this form we can specify a plasma flow via the Alfven 
Mach number Ma{A). Because A is constant on magnetic 
field lines this is also true for Ma and a. Physically this 
means that we can specify on which field lines plasma is 
flowing with a certain Mach number. It is in particular pos- 
sible to calculate separatrix field lines in the static case and 
specify plasma flow only on one side of this separatrix, e.g. 
to model plasma flow around a static magnetosphere or hel- 
met streamer configuration. 



Equation (22 1 is also equivalent to 



A = ±y" y/l-MA{aYda. 



(23) 



Some care has to be taken for multi-valued functions Ma [A) 
or A/4 (a), where one has to distinguish between the dif- 
ferent branches of solutions. This is, however, not a major 
problem and similar to the problem of multi-valued func- 
tions n(A) in the static Grad-Shafranov theory, which has 
been addressed in |WiegeImann et al[] ( |1998| l to model triple 
coronal helmet streamer configurations. 

2.1.1 Influence on the electric current density 

Inserting the ansatz B = Va x into (|4]i and by computing 
the Laplacian of a with equation ( [22| we find the connection 
between the Alfven Mach number and current density 



- MoJ^ = Aa 

= ±^^iVAf± 



-.AA. 



(24) 



The ± sign again reflects the freedom of the transformation 



relations Eqs. (22 1 and (23 1 with respect to the direction of 



the magnetic field and the symmetry of the Lorentz force. 
The first term in equation ( [24| corresponds to a current in- 
duced by the plasma flow and the second part modifies (en- 
hances) the static equilibrium current — Mo jz, static = — A^. 

For static potential fields this part of the current vanishes 
also in the stationary state with flow, i.e. as AA = only 



the first term on the right side of Eq. (24i contributes to the 



electric current density. We expand equation ( 24 1 for small 
Alfven Mach numbers Ma <C 1 which leads to 

-^,oJz^±iMA + ^Ml)M'A{VAf±{^ + ^Ml)AA,(25) 

and if we neglect all quadratic and higher terms in Ma, we 
find: 



^fioj, = ±MAM'^{VAf±AA. 



(26) 
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Inserting the definitions of the magnetic field and equilibrium 
current and using the definition il/^ V A = V A = VMa 



(this can be done also already in ( 24 1 and does not depend on 
the assumption of small Mach numbers), we get from ( [26] l 

-/ioj, =±(MAVAfA- VA + j2,static). (27) 

Consequently for small Alfven Mach numbers the equilib- 
rium currents are basically unmodified by the plasma flow 
and the induced currents depend linearly on the magnetic 
field strength, the Alfven Mach number and the gradient of 
the Alfven Mach number. For further approximations on the 
relative strength of the equilibrium currents and the flow in- 
duced currents we assume that the magnetic field equilibrium 
does change on a length scale /static and the plasma flow on a 
scale /flow, which allows us to roughly approximate the gra- 
dient and Laplacian: 



B 



(28) 



So to compare the relative strength of the two contributions 
we have to compare and and get as the ratio of in- 
duced and equilibrium current ■ Consequently we 
get (for slow plasma flows with Ma <C 1) only a significantly 
large induced current if the plasma flow changes on a much 
smaller length scale as the typical scale of the configuration 
/flow ^ /static- Such a situation is typically fulfilled at bound- 
ary layers, e.g., the magnetopause or the separatrix between 
open and closed field lines in coronal helmet streamers. 



3 Application to coronal helmet streamers and plas- 
moids 

In the following we provide some example solutions for 
MHD-equilibria with plasma flow. We construct these con- 
figurations by first solving the MHS problem pT| ) and then 
by transforming the resulting static flux function A{x,y) into 
the solution of the stationary problem a{x,y) with the help 



of equation ( 22 1. We prescribe the Alfven Mach number as a 



function of A in the useful form 

Ma (A) = Ml + (M2 - Ml ) tanh {d{A~Ac)), (29) 

where Mi, M2, d, and Ac are free parameters the scale 
on which the flow changes (doc l//flow is an inverse length), 
and the value of the separatrix field line Ac- The functional 
form of Ma in p9| ) has been chosen in order to provide the 
strongest flow gradient at the separatrix field line A^. Figure 
[T]top panel shows Ma as a function of A for Mi =0.4, M2 = 
0.8, d = 2, and Ac = 3. 

As an example we apply the transformation ( [22) l with Ma 
in the form of ( [29| to a homogeneous potential magnetic field 
B = BoBa: with Bq = 1. This is a simple 1-D equilibrium 



'Please note that the Alfven Mach number can become negative 
for plasma flows antiparallel to the field lines. 



with the static flux function A — Boy and all quantities (both 
in the static and stationary case) are only a function of y. Fig- 
ure [TJbottom panel shows the Alfven Mach number 
(dashed line), the streaming vector w{y) (dotted line) and the 
corresponding formation of a current sheet jz{y) at the sepa- 
ratrix field line Ac = 3 (solid line). 

In the following we study more sophisticated static equi- 
libria and their transformation to stationary incompressible 
MHD-equilibria. For a better visualization we will present 
examples of MHS equilibria that do not show the extremely 
small scale flows compared to the equilibrium current scales. 
We therefore have to use larger Alfven Mach numbers. 

3.1 Linear MHS equilibria 

To derive 2D static equilibria we solve the Grad-Shafranov 



equation (21 1 for a linear current. Such configurations have 
been studied for triple coronal helmet streamer configura- 
Wiegelmann ( 1998| l, but here we limit our research 



tions in 



to single helmet streamers and concentrate on the effect of 
plasma flow on open field lines. A linear current means that 



the function n(^) in equation (21 1 has the form Il{A) = 

2 ^ 

^A. In this case the Grad-Shafranov equation reduces to 
a linear Helmholtz equation 



(30) 



and can be solved by separation of variables. Let us remark 
that in general the static Grad-Shafranov equation has the 
form Il{A) = p{A) + -if, with the plasma pressure p and 
a magnetic shear field Bz in the invariant direction. In the 
case of p{A) = one obtains linear force-free configurations 
and else static equilibria. In both cases the electric current 
jz {A) = is linear in A. The particular choice c — cor- 
responds to current-free potential fields. As solution of (30i 
one gets by separation of variables 

A[x,y) = BqCxpI^ L^/*^°^\~Z^] < ' (31) 

A{x,y) — Bocos(^^^^^^ cos (^^^^j ioTc>k, (32) 



with J/ = ^/k'^ —c^ and oj = \/c 



Linear combinations 



of these particular solutions are also solutions of the linear 
Helmholtz-equation (30 1. These solutions were also stud- 



ied by Hood and Anzer ( 1990) modeling prominence arcades 
and [Petriej ( |2006) l modeling coronal loops. Here we consider 
only three particular cases with Bq — k — L — 1 and differ- 
ent values of c. The top panels of figures |2] |3] and |4] show 
magnetic field lines (equi-contour plots of the flux function 
A{x,y) for c — 0, c — 0.9, and c ~ 1.2, respectively). The 
case c = in figure |2] corresponds to a current-free poten- 
tial field. Introducing a moderate linear current with c< k 
leads to a stretching of the configuration (top panel in figure 
[3]) and a smooth electric current density distribution (second 
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panel in figure [3]l. For the case c> k, as shown in the top 
panel of figure 14] the magnetic topology changes and wet 
get plasmoid-like configurations, which, however, also have 
a smooth current density distribution in equilibrium (second 
panel). 

3.1.1 Transformation to stationary states 

We transform these static equilibria into stationary ones us- 
ing equation (22i and Af^ in the form (29 1 with d = 5, 
Ml = 0.0, Ma = 0.8, A^. = 0. We choose A/a > Mi, in or- 
der to prescribe a plasma flow on the outside (x — ±0.5) of 
the configuration, where the flux function becomes negative. 
Inside the streamer and plasmoid (—0.5 < x < 0.5) the flux 
function is positive. The plasma flow is chosen in order to 
be maximal at the open separatrix field fine and the plasma 
is basically at rest inside the configuration, where the mag- 
netic field lines are closed (bottom panels in figures [2j [3] |4|i. 
With the current transformation equation (|24]) we compute 
the total current density: In the case of the static potential 



field currents only occur from the first term in ( 24 1 and are 



driven by the gradient of the plasma flow (center panel in fig- 
ure |2]i. As a consequence of the sharp gradient in the flow, a 
thin current sheet forms at the separatrix (x = ±0.5). Such a 
thin current sheet forms also additional to the smooth equi- 
librium currents for the linear current cases (third panels in 
figures [3] and |4]). In regions with weak or no plasma flow 
the equilibrium current does basically not change, whereas 
a current sheet forms at the separatrix field line. Due to the 
strong gradient in the plasma flow this induced current sheet 
is much thinner and the current density is higher than in the 
equilibrium current. The equilibrium currents are strongest 
in the center of the configuration. 



4 Application to magnetosplieres 

4.1 Non-linear Grad-Shafranov equation, Liouville 
equation 

In the static case without magnetic shear field the Grad- 



Shafranov-equation (21 1 reduces to 



(33) 



where p{A) is the plasma pressure. Under the assumption 
of a local thermodynamical equilibrium the plasma pressure 
function can be derived from kinetic theory in the form (see 
|Schindler| [20061 for details) 



piA) = -pexp{~2A/A), 



(34) 



where A and p are normalization constants. This leads to an 
equation in the form 



AA = Aexp(-cA), 



with constants 

P 



A 



Mo 



A 



A' 



(36) 
(37) 



and the typical lengthscale 1, and the typical magnetic field 
B, defined by 



Ac 



2p 



2 



(38) 



(39) 



I 

Before we continue to compute analytical and exact 2D so- 
lutions, we first present a well-known case in ID, namely the 
Harris-sheet, to explain how the transformation from MHS 
to incompressible SMHD works. 

4.2 ID Harris-sheet 

A well known, ID equilibrium current sheet solution of the 



Liouville's equation (35i is the Harris-sheet (Harris 19621. 



A Harris -sheet like force-free equilibrium has recently been 



found by Harrison and Neukirch| (2009) . The Harris sheet 
is a ID-solution, where all quantities depend only on the y 
coordinate, and is given by 



A{y) ^ A\ncosh{y /I) . 



(40) 



The static equilibrium quantities for the Harris-sheet are 
shown with solid lines in figure [5] The top panel shows the 
flux-function A{y), the second panel the corresponding mag- 
netic field Bx and the bottom panel the equilibrium electric 
current density. 

4.2.1 Transformation to stationary states 



We use the transformation equation ( 22 1 with the Mach num- 
ber profile (29 1 and Af„iax = 0.5, Mmin — 0, Ac — 2, and 



d = 5 to derive a stationary equilibrium with plasma flow. 
The stationary solutions are shown with dashed lines in fig- 
ure|5] The top panel shows the flux function a{y), the second 
panel the magnetic field Br^ and the bottom panel the electric 



current density (y) as computed with equation (24 1. Addi 



tionally we present the current approximation for Ma ^ 1, 



(35) 



as computed with equation ( 26 1 with dotted lines in the bot- 
tom panel of figure |5] By comparing the static (solid lines) 
and stationary (dashed) quantities, one can see that both 
quantities only differ in the region were the flow gradient is 
high (see third panel). In these regions, where the static equi- 
librium is separated from the stationary flowing plasma in a 
thin layer, current sheets form. The spatial scale of these lay- 
ers is significantly smaller than the typical length scale of the 
equilibrium current. The electric current density approxima- 
tion for small Mach number (dotted in bottom panel) shows 
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reasonable agreement with the exact (dashed) solution, even 
for the not very small maximum Mach number Ma —0-5 
used here. 

4.3 Exact 2D magnetospheric equilibrium 



As sho wn by |LiouyiUe| ( [T853] l, |Bandle| ( [T975] l, |Birn et al.| 
( 1978 1 equation |35| can be written as: 



dudu 



■ Aexp(— cA), 



(41) 



with u — x + iy and u — x — iy, and this equation has the 
general solution 



A(u,u) — - In 
c 



2, 1 + f 



,0*1 

I du I 



= A\n- 



1+- 

4 








a* 

9it 





(42) 



Every analytic function '^{u) generates a solution of (35 i. 
The Liouville equation has also its applications outside 
plasma physics and solutions in the form of ( [42| have been 
used, e.g., by ,Schmid-Burgk| ( |1967| ) to investigate a self- 
gravitating gas layer. 



Schindler and Birn ( 2004 1 found magnetospheric solutions 



with the Ansatz 




(43) 



which leads to the solution class of ( |42| i of the form (see 
[Schindler and Bim] |2004| for details and discussion of the 
static equilibrium). 



A(x,y)/i = ln 



cosh 







(44) 



with r = y/x^~+y^, where the coordinates here are nor- 
malized on /. For e — !• oo the Ansatz for ^I^ produces the 
Harris-sheet solution. The top panel of figure |6] shows the 
corresponding magnetic field lines as equi-contour plots of 
A{x,y) for e = 1. For the transformation we used the param- 
eters Ml = 0.0, M2 = 0.95 and d = 2 for the Mach number 
profile in Eq.(|29]). The second panel contains the equilibrium 
current density —jz{x,y)- 

4.3.1 Transformation to stationary states 



We use the transformation equation ( 22 1 with the Mach num 

'mm — 0, Ac = 



ber profile (|29]l and A/,„ax = 0.95, Afmin = 0, Ac = 0, and 
d = 2 to derive a stationary equilibrium with plasma flow. 
The chosen profile for the flow is smoother, as in previous 
examples and the maximum Mach number is higher. As 
a consequence we observe two additional current sheets in 



the stationary current distribution as computed with ( 26 1 and 
shown in the third panel of figure|6] The thickness of the two 
current sheets are located in the region with the plasma flow 
gradient, shown in the bottom panel. The smoother profile of 
Ma{A) results also in smoother induced current sheets. 



5 Conclusions 

Within this work we studied the relation between plasma 
flow gradients and current sheets in space plasma. To high- 
light the influence of stationary flows on static MHD equi- 
libria, we neglected compressibility effects and used the as- 
sumption of field-aligned, incompressible stationary flows. 
These assumptions imply an analogy between magnetic field 
and velocity field as well as an analogy between MHS and 
incompressible SMHD: The assumption of incompressibility 
allows us to transform magnetostatic equilibria into station- 
ary ones by using a non-canonical transformation. We find 
that the occurrence of flow driven current sheets is closely 
related to the gradient of the plasma flow or, to be precise, to 
the gradient of the Alfven Mach number perpendicular to the 
magnetic field lines. Along the field lines the Mach number 
is always constant for incompressible stationary flows. As 
the gradients in the Alfven Mach number can be very large, 
because the typical length scale of the flow is smaller than 
the length scale of the magnetic field, the occurence of cur- 
rent sheets is correlated with the appearance of vortex sheets. 
Such configurations can be closely connected to local break- 
downs of the frozen-in-flux theorem, or shortly to magnetic 
reconnection ( |Eyink and Aluie[|2006| l. Magnetic reconnec- 
tion in turn plays a major role in eruptive space plasma pro- 
cesses like magnetospheric substorms or solar flares. 

In principle it is possible to compare our theoretical in- 
vestigations on the relation between plasma flows and cur- 
rent sheets with observations, in particular in the magneto- 
sphere where in-situ measurements are available. One pos- 
sibility is using magnetic field and particle data from the 
CLUSTER-mission, which are carried out simultaneously 
with four spacecraft. Such multi-spacecraft measurements 
(with distances between the spacecraft in the range of about 
50-10000 km) of the magnetic field allow also the estima- 
tion of electric currents. By taking moments of the particle 
data it is possible to compute plasma quantities like density, 
pressure and the plasma flow velocity. These combined mea- 
surements allow at least to estimate gradients in the plasma 
flow and the thickness of current sheets. A limitation is that 
structures smaller than the distance of the Cluster-spacecraft 
cannot be spatially resolved, which implies that flow gradi- 
ents could be steeper and the current sheets thinner as com- 
puted from the measurements. A comparison of data with 
our model, which relates flows and flow gradients to current 
sheets, will allow to investigate how consistent different ar- 
eas in the magnetosphere can be described under the assump- 
tion of stationary incompressible MHD. 
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Fig. 2. 2D potential field configurations in the form A{x,y) = 
-Boexp(— fc7ri//L)cos(— x/fc'^ — c^ttx/L), with Bq — k — L = 1, 
c — 0.0 and the Alfven Mach number profile function \29\ with 
d = 5, Ml = 0.4, M2 = 0.0, Ac = 0. This profile has a steep gradi- 
ent at the boundary between open and closed field lines. Top panel: 
magnetic field lines (contour lines of A{x,y)), center: formation of 
flow driven current sheets Jz{x,y), Bottom: Alfven Mach number 
MA{x,y). 




Fig. 3. Same transformation as in figure|2] including a linear MHS 
current Jz{A) = c - A, c = 0.9. Second and third panel show the 
split of linear (MHS, without flow) and the full current (with flow). 



Fig. 4. Same transformation as in figure |2] but including a linear 
MHS current in the form Jz^A) = c- A with c= 1.2. Here c> k. 
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Fig. 5. From top to bottom we show the flux-functions A(y),a{jj), 
the magnetic field Bx{y), the Alfven Mach number MA{y) and the 
electric current density jz{y)- Solid lines correspond to the mag- 
netostatic case and dashed lines to stationary MHD with the profile 
( (29} and Mi = 0, M2 = 0.5, Ac = 2 and d = 5. In the bottom panel 
we show additional dotted the approximation for Ma ^ 1 as com- 
puted with \27\ . 
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Fig. 6. Same transformation formula as figure|2] but for a nonlinear 
current in the form Jz{A) oc exp{—cA). For the transformation we 
used A/i = 0.0, M2 — 0.95 and d = 2. For a better visualization we 
show the negative current — jz {x,y) in the second and third panel. 



